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TWO-DIMENSIONAL COMPRESSIBLE FLOW IN CENTRIFUGAL 
COMPRESSORS WITH STRAIGHT BLADES 

By John' D. Stanitz and Gaylord O. Ellis 


SUMMARY 

Six numerical examples are presented for steady, two- 
dimensional, compressible, nonviscous flow in centrifugal 
compressors with thin straight blades, the center lines of which 
generate the surface of a right circular cane when rotated about 
the axis of the compressor. A seventh example is presented for 
incompressible flow. The solutions were obtained in a region of 
the compressors, including the impeller tip, that was considered 
to be unaffected by diffuser vanes or by the impeller-inlet configu- 
ration. Each solution applies to radial- and mixed-flow com- 
pressors with various cone angles but with the same angle 
between blades on the conic flow surface. ( The solutions also 
apply to radial- and mixed-flow turbines with the rotation and 
the flow direction reversed.) The effects of variations in 
the following parameters were investigated: ( 1 ) flow rate, 
iff) impeller-tip speed, ( 8 ) variation of passage height with 
radius, and (f) angle between blades on conic flow surface. 
The numerical results are presented in plots of the streamlines 
and constant Mach number lines. 

Correlation equations are developed whereby the flow condi- 
tions in any impeller with straight blades can be determined 
(in the region investigated by this analysis) for all operating 
conditions. As examples of the information provided by the 
correlation equations, the velocities along the blade surfaces are 
presented for a wide range of impeller-tip Mach number, flow 
coefficient (flow rate), and angle between blades on the conic 
flow surface. 

INTRODUCTION 

At the present time, the design of centrifugal compressors 
is primarily an art rather than a science. Little detailed 
knowledge of flow conditions within the compressor on which 
to base a rational design is available. If these flow condi- 
tions could be determined, design methods might be devel- 
oped for centrifugal compressors with higher aerodynamic 
efficiency and better over-all performance. For example, 
the compressor efficiency would be improved if favorable 
velocity distributions (with respect to boundary-layer growth 
and separation) could be obtained along the flow surfaces by 
proper design of the compressor. 

For a given set of design and operating parameters, the 
velocities and pressures within the compressor depend on the 
three-dimensional flow path and on the fluid properties 
(compressiblity and viscosity). A complete analysis of the 
flow must include all these factors. If flow conditions are 


essentially uniform in one direction, however, the flow is 
adequately represented by a two-dimensional analysis in 
which the fluid is considered inviscid but compressible. 
Viscosity of the fluid is unimportant except within the bound- 
ary layer along the flow surfaces, and this boundary layer 
is thin provided favorable velocity distributions exist 
within the compressor. On the other hand, compressibility 
of the fluid is important in centrifugal compressors because 
the large pressure ratios per stage result in density changes 
that affect fluid velocities, streamlines, and so forth. 

In a previous report (reference 1), a general method of 
analysis is developed for stead\ r , two-dimensional, compres- 
sible flow through radial- and mixed-flow compressors and 
turbines in which the center line of the passage generates 
the surface of a right circular cone when rotated about the 
axis of the machine. The two-dimensional flow pattern is 
considered to lie upon the surface of this cone. 

In the present report, these analytical methods are applied 
to investigate the flow conditions within a certain region of 
radial- or mixed-flow compressors (and turbines) with thin 
straight blades lying on conic radii (elements). The region 
investigated includes the impeller tip and is that region 
considered to be unaffected bj' the inlet configuration of the 
impeller and by the diffuser vanes; that is, the impeller inlet 
and the diffuser vanes (if any) must be far enough removed 
from the region investigated not to affect the flow appre- 
ciably in that region. Straight blades were selected because 
they were considered the most representative blade shape 
now in use for aircraft centrifugal compressors. 

The purpose of this analysis was to determine the effect 
of operating and design variables (impeller-tip speed, com- 
pressor flow rate, variation of passage height with radius, 
and angle between blades on conic flow surface) on flow 
conditions within the region investigated. From this 
information, limitations can be placed on the operating and 
design variables if certain flow conditions are desired within 
these regions. For example, if from boundary-layer con- 
siderations maximum rates of deceleration of the relative 
velocity are specified along the flow surfaces, the results of 
this analysis can be used to determine limiting values of 
impeller-tip speed, compressor flow rate, variation of passage 
height with radius, and angle between blades on conic 
flow surface. 

The theoretical investigation presented herein was con- 
ducted at the NACA Lewis laboratory during 1948. 
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METHOD OF" SOLUTION 
EQUATIONS 

A general analysis is developed in reference. 1 for steady, 
two-dimensional, compressible flow in. compressors and 
turbines with arbitrary blade shapes and arbitrary varia tions 
in the passage height. The analysis isTimited to radial- and 
mixed-flow compressors and turbines in. which the center 
line of the passage generates the surface of a right circular 
cone when rotated about the axis of the machine (fig. 1). 



Figure 1.— Fluid particle on rotating coordinate system ot impeller. Center line of flow 
passage generates surface, of right circular cone with cone angle a. 


The two-dimensional flow pattern is considered to .lie on the 
surface of this cone. The method of analysis is applied herein 
to radial- or mixed-flow compressors with thin straight 
blades lying on conic radii (elements). The equations de- 
veloped in reference 1 are presented in this section with a 
brief discussion of the coordinate system and the assumptions 
and the limitations of the analysis. 

Coordinate system. — A developed view of the conic flow 
surface generated by the passage center line (fig. 1) is shown 
in figure 2. The dimensionless conic„coordinates of a fluid 
particle on the conic flow surface are R and 9. (All symbols 
are defined in appendix A.) The coordinate d is considered 



Figure 2 .— Fluid particle op developed view ot conie surface. Dimensionless coordinates 
relative to Impeller R, 6, and H\ tangential and radial compononts of velocity relative to 
impeller u and v, respectively. 


positive in the counterclockwise direction. The conic- 
radius ratio R is defined by 

where r is the conic radius (distance along conic element 
from apex of cone) and the subscript T refers to the impeller 
tip. The coordinate system (R,d) rotates in the positive 
direction., of 6 with the angular velocity u of the impeller. 
The passage-height ratio H, in the direction normal to jhc 
conic flow surface (fig. 2), is a continuous function of the 
conic-radius ratio R. 

H =r= m - - (2) 

where h is the passage height at any conic-radius ratio R. 

Assumptions and limitations. — This analysis assumes that 
flow conditions are uniform across the passage normaT to 
the conic flow surface; that is, the flow varies only along the 
conic flow surface. In order to satisfy this assumption, it 
is necessary that: (1) the gradient of h with respect to r be 
small; and (2) the cone angle a be sufficiently close to 180°. 
The allowable variation in a from 180° will depend on the 
ratio h/r and on the desired accuracy. For the hypothetical 
limiting case in which h[r approaches zero everywhere along 
the conic flow surface, the analysis is accurate for all values 
of or. _ • 

Velocity-ratio components. — The fluid particle on the 
developed conic flow surface in figure 2 has a relative 
tangential-velocity ratio U and a radial- (along conic element) ' 
velocity ratio V. These velocity ratios are defined by 

... . . .. (3a) 

and 

V~ (3b) : 

c o 



TWO-DIMENSIONAL COMPRESSIBLE PLOW IN CENTEIFU GAL COMPRESSORS WITH STRAIGHT BLADES 


where 

u tangential component of velocity relative to impeller 
in positive direction of & (fig. 2 ) 
v radial (along conic element) component of velocity 
(fig- 2) 

c local speed of sound 
Subscript: 

o absolute stagnation condition in region of uniform flow 
upstream of impeller 

Stream function f . — A dimensionless stream function \p 
satisfies the continuity equation if defined as 

§H VHE (4a) 

and 

« b > 

where p is the weight density. 

The stream function \p is constant along the blade surface. 
If \p and 6 are assigned values of zero along the positive 
blade surface (the blade surface in the positive direction 
of 0), the value of \p along the negative blade surface (the 
blade surface in the negative direction of 0 ) is given by 


^„=ipA6 


( 5 ) 


where A0, the angle between blades on the conic flow surface, 
is given by 

, „ 2 ~ . a 

A0=-g sin y (6) 

and where 


W 

PqCLtCo 


( 7 ) 


' where 


a j* — B (A0)r fK j. 


(7a) 


where 

<p compressor flow coefficient 

IF total compressor flow rate 

a annular flow area (normal to conic flow surface) 

B number of passages (or blades) 

Subscript: 

n negative blade surface (blade surface in negative 
direction of 6) 

Differential equation. — The differential equation for the 
stream-function distribution in compressors with thin 
straight blades lying on conic radii (elements) is obtained 
from equation (14) of reference 1 in terms of transformed 
coordinates £ and 17 . 


2 M, 


t- m* 

Po 


exp [(m- 


-2)^0]-^+^ ^ 


<AP 

£>4 



m(A0) || 


d 



( 8 ) 
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where the impeller-tip Mach number M r is defined as 


M t 


m r T sin — 

~C~o 


( 9 ) 


The transformed coordinates { and ij, which for thin straight 
blades lying on conic radii are related to R and 6 by 


t log, R 

A 9 


( 10 a) 


0 


V== A0 


( 10 b) 


have been introduced because they result in parallel blades 
in the transformed plane. Such a transformation is desirable 
because it simplifies the solution of the differential equation 
by relaxation methods. 

In equation ( 8 ), the passage-height ratio H is assumed to 
vary with the conic-radius ratio according to 

H=R m ( 11 ) 

where m is an arbitrary exponent. (For m=0, the blade 
height remains constant and for m — — 1 . 0 , the flow area normal 
to the conic flow surface remains constant.) 

In order to solve equation ( 8 ), it is necessary to know the 
density ratio, which is related to the impeller-tip Mach 
number and the relative velocity ratio Q by (equation ( 11 ), 
reference 1 ) 


1 

(i2) 

where y is the ratio of specific heats and 

(13) 


and the absolute whirl of the fluid upstream of the impeller is 
assumed to be zero. The velocity Q (multiplied by p/p„) is 
in turn given by equations (4), (10), (11), and (13) as 


p 0 v (AS) exp [(m+ 1)$A0] L\c>f / ^V>V ) J 


(14) 


Equations ( 8 ), ( 12 ), and (14) provide three equations with 
three unknowns p/p„, Q, and 1 p. 

NUMERICAL PROCEDURE 

The system of equations ( 8 ), (12), and (14) is solved by 
relaxation methods to obtain the stream-function distribution 
within the compressor. From this distribution, the velocity 
components and other conditions can be determined using 
equations (4) and so forth. Detailed outlines for the 
numerical procedure are given in references 1 and 2 . The 
procedure is briefly sketched herein. 
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Outline of procedure, — In order to solve the system of 
equations by relaxation methods, the following procedure is 
employed: ”... 

(1) Equations (8) and (14) are changed to finite-difference 
form (to be discussed). 

(2) Values of 4 are specified on the boundaries of the flow 
region. 

(3) Values of \p are estimated at equally spaced points of a 
grid system within the boundaries of the flow region. 

(4) The preceding estimated values of 4/ are adjusted 
(relaxed) by the relaxation process until they satisfy equa- 
tion (8) in finite-difference form. 

(5) The boundary values of ^ in the vaneless diffuser are 
adjusted to satisfy the Joukowski condition for tangency of 
flow leaving the tip of the impeller blade. 

(6) After the Joukowski condition has been satisfied, the 
grid spacing is reduced near the impeller tip in order to obtain 
detailed knowledge of the flow characteristics in this region 
where conditions are rapidly changing. 

Finite-difference equations. — Equation (8) is changed to 
the following finite-difference form (reference 1) : 

^i+^+^+»4-4^- ( ^7^ (log«g-log.g)-- 

2M t £(A 6) 2 exp [(m+2)£A0] b 2 = R (15) 

Po ■ ......... 

where ' 

b grid spacing (fig. 3) , 

R residual (error) due to estimated values of 4 / used dur- 
ing relaxation solution.... 



Figube 3.— Sample grid shoving grid spacing b and numerical subscript, convention for 
adjacent grid points. 


Subscripts: 

1, 2, 3, and 4 quantities at adjacent grid points as defined 
V in figure 3 (Quantities without numerical 

subscripts refer to grid point at which R is 
being computed.) 

From equation (12), the natural logarithm of the density 
ratio (required for the solution of equation (15)) is plotted 
as a function of the flow-rate ratio Qp/p 0 in figure 4. The 
flow-rate ratio is obtained from equation (14), which in 
finite-difference form becomes 

P Q . Wi— fo)*+(fo— ifrO 2 ]* (16) 

' p 0 ' 2b {AS) exp [(m-j-l)£A0] 



Figure 4. — Natural logarithm o£ density ratio log. pip. as function o£ flow-rato ratio Qplp- lot 
. various values of RMt. Equation (12); ratio of speeifio heats y, 1.4. 
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Fiqure 5.— Relaxation grid in plana o£ transformed coordinates. 
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Boundary values of ip . — The stream function ip is constant 
along the blade surfaces. If the value of p is arbitrarily 
set equal to zero along the positive blade surface, the value 
along the negative blade surface is given by equation (5). 
Values of ip along the left boundary (fig. 5) are obtained 


Flow region Investigated 



from the simplified analysis of reference 1 , which assumes 
that the flow is parallel to the blade surfaces. Values of 
ip along the boundaries of constant y in the vaneless dififuser 
are estimated by a one-dimensional analysis of the flow 
(reference 1 ), which assumes the moment of momentum to 
be constant. (Mote that these boundary values of ip in the 
diffuser are only approximate and must therefore be relaxed 
along with the values of ip at interior grid points, as indicated 
in reference 1 .) At any value of £ in the vaneless diffuser, 
the values of p along the boundary y = 1.0 are p„ greater 
than the values of ip along the boundary y=0. The values 
of ip along the downstream boundary of the diffuser (fig. 5) 
vary linearly. 

Estimated interior values of ip . — Correlation equations 
developed in this report (and presented in a later section) 
provide a first approximation for \p at grid points within the 
impeller passage. The error in these estimated values of 
p will, in most cases, be less than ±1 percent of 

Relaxation process.' — The residuals R that result from the 
estimated interior values of 1 p are computed at each grid 
point by equation (15). These residuals are then reduced 
(relaxed) by suitable changes in the values of ip. The detailed 
procedure is given in references 1 and 2 . 

Joukowski condition. — The Joukowski condition requires 
that the streamline leaving the blade surface be tangent to 
the blade tip. This condition generally is not satisfied by che 
initial relaxation solution because for this solution the bound- 
ary values of ip along the downstream boundary in the vane- 
less diffuser (fig. 5) are obtained from the estimated variation 
in ip along the boundaries of constant 13 in the diffuser (dis- 
cussed in the section Boundary values of if). In order 
to satisfy the Joukowski condition, the values of p along 
the downstream boundary must all be changed by the same 
required amount and this change in the boundary values of \p 
results in changes in ip at each of the interior grid points. 
The manner in which these changes in ip are made to satisfy 
the Joukowski condition is discussed in reference 1. 


RESULTS 

Seven numerical examples are presented. One of these 
examples has been selected as the “standard” and in each 
of the remaining examples one parameter is varied from the 
standard conditions as shown in the following table: 


Example 

V 

Mr 

m 

Afl 

(deg) 

Type of flow 

Standard. 

0.5 

L5 

- 1.0 

12 

Compressible ( 7 = 1 . 1 ) 

1 

.7 

L5 

-1.0 

12 

Compressible (7=L4) 

2... 

.9 

L5 

-1.0 

12 

Compressible ( 7 = 1 . 4 ) 

3 

.5 

t.o 

-LG 

12 

Compressible ( 7 -1.4) 

4 

.5 

L5 

- 1.1 

12 

Compressible (7=L4) 

5 

.5 

1.5 

-1.0 

IS 

Compressible ( 7 = 1 . 4 ) 

6 

.5 

1.5 

—LG 

IS 

Incompressible 


By comparing the solution for the standard example with 
one of the nonstandard examples, the effect of the change 
in a single design or operating parameter on flow conditions 
in the compressor can be determined. 

These examples are for impellers having straight blades 
(fig. 6 ). The solutions were obtained in a region of the 
compressors (including the impeller tip, fig. 6 ) that was 
considered to be unaffected by the inlet configuration of the 
impeller and by the diffuser vanes; that is, the diffuser vanes, 
if any, must be far enough removed from the impeller not to 
affect the flow region being investigated. Each solution 
applies, within the limitations imposed by the assumption of 
two-dimensional flow, to radial- and mixed-flow compressors 
(and turbines) with conic flow surfaces having various cone 
angles a but the same angle between blades on the coni c flow 
surface A0 (reference 1). 

The numerical results are presented in plots of the stream- 
lines and constant Mach number lines. 

Streamlines. — The streamline configurations (relative to 
the impeller) for the seven examples are shown in figure 7. 
The streamlines are designated in such a manner (P/p n ) that 
the value of a streamline indicates the ratio of the flow that 
lies between the streamline and the positive blade surface to 
the total flow in the passage. For a given density ratio, the 
streamline spacing is indicative of the velocities relative to 
the impeller, with close spacing indicating high velocities and 
wide spacing indicating low velocities. 

The streamlines for the standard example are given in 
figure 7 (a). For the design and operating conditions of this 
example, an eddy has begun to form on the positive blade 
surface. This eddy results from negative velocities on and 
near the blade surface. The fluid in this eddy rotates 
(relative to the impeller) in the opposite direction to that of 
the impeller so that the absolute motion of the fluid is irrota- 
tional. In actual practice this eddy is probably unstable 
and it is desirable to eliminate the eddy by proper changes in 
the design and operating conditions of the compressor. 
From an inspection of figure 7, it appears that the eddy can 
be reduced or eliminated by increasing the flow coefficient <p 
(figs. 7 (b) and 7 (e)), decreasing the impeller-tip Mach 
number M r (fig. 7 (d)), decreasing the angle between 
blades on the conic flow surface Ad (fig. 7 (f)), and usings 
incompressible fluids (fig. 7 (g)). The eddy in figure 7 (d) 
(impeller-tip Mach number of 2.0) is especially interesting 
because it occupies more than one-half of the available flow 
area at a radius ratio of 0.90. 




(a) Standard example: flow coefficient OJS; 
impeller-tip Mach number Mr, 1.5; 
constant flow area (ro,— 1.0); angle be- 
tween blades on conic flow surface A9, 12°; 
compressible flow Ct. 1.4). 


(b) Example 1 : flow coefficient *>» 0.7; other 
parameters same ns standard example. 



Conic-radius ratio, 



(c) Example 2: flow coefficient tf, 0.9; other (d) Example 3: impeUer-tip Mach number 

parameters same as standard example. Mr, 2.0; other parameters same as stand- 

ard example. 

Fxgube Continued. Relative streamlines’for flow throngh'centrifugal compressor with 
straight blades. Streamline designation indicates ratio of flow that lies between streamline 
and positive blade surface (right side of passage) to total Sow in passage. 
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u 


Relative 

streamlines 



(e) Example 4: varying flow area (m,— 1.4); 
other parameters same as. standard ex- 
ample. - 

Figuee 7.— Continued. Relative streamlines for flow through centrifugal compressor with 
straight blades. Streamline designation indicates ratio of flow that lies between stream- . . 
line and positive blade surface (right side of passage) to total flow in passage. 

The exponent m determines the change in flow area 
through the impeller and the diffuser. For m equal to —1.4 
(example 4), the flow-area ratio (flow area divided by flow 
area at impeller tip) at a radius ratio of 0.5 is 1.32 compared 


to a ratio of 1.0 for the sta ndard example (in which in 
equals —1,0). The effect of m on the streamline configura- 
tion in the flow region investigated is inappreciable (compare 
figs. 7 (a) and 7 (e)). The. reason for this small effect is that 
the flow areas at the impeller tip are the same (because <p is the 
same) and for reasonable values of m are not much different 
anywhere in the flow region investigated (which is in the 
vicinity of the tip, fig. 6). If the area ratio 1.32 had been 
obtained by maintaining equal flow areas for both impellers 
at the radius ratio 0.5 within the impellers and decreasing 
the area at the impeller tip in one case, the streamline con- 
figuration would be greatly affected. This effect, however, 
would not he the result of the, change in m but rather of the 
change in \p (which results from the change in a r ). 

The large effect of compressibility on the streamline 
configuration (and therefore on the other flow conditions) is 
shown by a comparison of figures 7 (f) and 7 (g) . The large 
eddy that exists for compressible flow completely disappears 
for incompressible flow. 

Mach number lines. — Lines of constant Mach number 
relative to the impeller are shown for the seven examples in 
figure 8. It should be noted that .the. Mach number in the 
incompressible solution (fig. 8 (g)) is a fictitious quantity 
that is equal to the fluid velocity q divided by a constant 
that is equal to the inlet stagnation speed of sound c 0 of 
whatever compressible solution with which the incompressi- 
ble solution is being compared. This inlet stagnation speed 
of sound is also contained in the definitions of <p and M T , so 
that for the incompressible solution qjc 0 , <p, and M T vary 
inversely with the assigned value of c 0 , but ratios of those 
parameters are unaffected. 

The standard example is given in figure 8 (a). The gen- 
eral characteristics of these plots are similar. The velocities 
(as indicated by the. Mach number lines) along the positive 
blade surface are low; the velocities along the negative blade 
surface are high; and the velocities become equal on the 
positive and negative blade surfaces at the blade tip (as 
required by the Joukowsld condition), The maximum 
Mach number occurs on the negative blade surface at a 
radius ratio, well within the impeller and the flow decelerates 
along the .blade surface from, this point to the blade tip. 
This deceleration, which for impellers with straight blades 
and the usual type of area variation with radius ratio becomes 
rapid near_the blade tip, is conducive to boundary-layer 
separation, which lowers the compressor efficiency. 

If the Boundary-layer wake in the vaneless diffuser is 
neglected, the velocities become essentially uniform at a 
radius ratio of approximately 1.10 for A6 equal to 12° 
(figs. 8 (a) to 8 (e)) and at a radius ratio of approximately 1.15 
for A 6 equal to 18° (figs. 8 (f) and 8 (g)). These radius 
ratios and their corresponding angles are equivalent to a 
value of £ approximately equal to 0.45. Flow conditions in 
the vaneless, portion of the diffuser immediately following 
the impeller therefore become essentially uniform at a value 
of £ approximately equal to 0.45. 
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<a 


Relative 
s t ream lines 



(f) Example 5: angle between blades on 
conic flow surface AS, 18°; other parameters 
same as standard example. 


Fig t : r s 7.— Concluded. Relative streamlines tor flow through centrifugal compressor with 

and positive blade surface (right side 

For the compressible-flow examples, the relative Mach 
number at the impeller tip is low (even for large values of <p) 
because of the high density ratios, which result from the 
high impeller-tip speed. These velocities would be con- 
siderably higher if the effective flow area were reduced by 
boundary-layer separation, which might be expected in 
practice. 

From an inspection of figures S (a) to 8 (g), it appears 
that the maximum Mach number (on the negative blade 
surface) is increased by increasing the flow coefficient <p 
(figs. 8 (b) and 8 (c)), is apparently not much affected by- 
increasing the impeller-tip Mach number M T (fig. 8 (d)) or 
by' changing the exponent m (fig. 8 (e)), and is increased 
by increasing the angle between blades on the conic flow 
surface A8 (fig. 8 (f)) or by changing to an incompressible 
fluid (fig. 8 (g)) (in which the Mach number is a fictitious 
quantity, as previously indicated). 


V 

Relative 



(g) Example 6: incompressible Quid; other parameters same as example 5 (fig* 7(0). Note 
that * and A/r are based upon constant c ff , which if considered equal to speed of sound at 
inlet conditions of example 5 enables comparison of compressible (example 5) and incom- 
pressible (example 6) solutions for same impeller-tip speed, weightrfiow rate, and so forth. 

tight blades. Streamline designation indicates ratio of flow that lies between streamline 
assage) to total Sow in passage. 

Slip factor. — The impeller slip factor is defined as the 
ratio of the average absolute tangential velocity of the fluid 
at the impeller tip to the tip speed of the impeller. The 
slip factor has been computed for each of the seven examples 
by methods given in reference 1. The resulting slip factors 
are given in the following table: 


Example 

Nonstandard 

parameter 

Slip 

factor 

Standard 


0.934 

1 

**=0.7 

.937 

2 

•-=0.9 

.908 

3 

Mr-2.0 

.935 

4 

i»=— 1.4 

.934 

5 

Afl = 18° 

.899 

6 

p/jb.-LO* 

.892 


‘Also, AS=18°. 


It appears that the only variable investigated that affects 
the computed slip factor is the angle between blades on the 
conic flow surface A8. In particular, it will be noted that 
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(a) Standard example: flaw coefficient <p, 0.5; 
impeller-tip Much number Mr, 1.6; 
constant flow area (in,— 1.0); angle be- 
tween blades on conic flow surface A3, 12°; 
compressible flow (v, 1.4), 
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Relative- 



(b) Example 1: flow coefficient V, 0.7; other 
parameters same as standard example. 


Figure 8.— Lines of constant Mach number relative to impeller. 
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(c) Example 2: flow coefficient 0.9; other (dl Example 3: impeller-tip Mach number 

parameters same as standard example. AJr. 2.0; other parameters same as stand- 

ard example. 

Figure 3. — Continued. Lines of constant Mach number relative to impeller. 
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ie) Example 4: mylngflow area (m,— 1.4); 

other parameters same as standard 

example. =--• * " — • ' --- - - - ---- ■- 

Figure 8.— Continued. Lines of qpnstant .Mach number relatiye to impeller. 


the slip factor is approximately the same for compressible and 
incompressible flow (compare examples 5 and 6), although 
the streamline configurations for the two examples are very 
different (compare figs. 7 (f) and 7 (g)). 


CORRELATION OF RESULTS 

Correlation equations are presented whereby the flow 
conditions (U, V, and so forth) within any impeller 
with straight blades can be determined (for the flow region 
investigated herein, see. fig. 6) from the flow conditions for 
the standard solution of this report. These correlation 
equations are developed in appendix B and the flow condi- 
tions forjhe standard solution are given in tables I to III. 

The correlation, equations are developed in terms of the 
transformed coordinates (£ and rj). The dimensionless conic 
coordinates R and 6 are related to these transformed co- 
ordinates by 

R=ezp(%A6) (B3) 

and _ 

' 8=--yAB (B4) 

Stream-function ratio >/#„.— The stream-function ratio \ 
varies across the impeller passage from 0 along the posi- 
tive blade surface of one blade to 1.0 along the negative 
blade surface of the next blade. At any given point (£ and ij) 
within the impeller, the value of the stream-function ratio 
for any impeller (with straight blades) and for any operating 
condition can be estimated by the following correlation 


equation (appendix B) : 

_ 

+ m 

1 

i_^_r 

r 


"VCn— 1) {y [A8-(A6),} + {R-R,UMrAd) 

,}) (B16) 

where the prime indicates the estimated value of flow condi- 
tion (stream-function ratio in this case) at a given point 
(£,17) and the subscript s indicates the standard value of 
flow condition at the same point. Also, 

- w j * 

. (M r AS), 

(B2) 

Pbrh 

Pa 

- (B9) 


and, if the .absolute whirl ahead of tho impeller is zero, 



where the subscript m indicates tho mean value at a given 
radius ratio R (that is, at a given value of £, which is related 
to R by equation (10a)). 

The estimated, values of the stream-function ratio \ /-'M. 
obtained from tho correlation equation (B16) are compared 
in figure 9 with the relaxation- values of \ pty* obtained for the 
numerical examples of this report. Values are plotted for 
every other grid point across the passage at- values of { indi- 
cated by the symbols. Perfect correlation corresponds to 
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Ui 


Ref a f/ve 
Mach number 
0.55 



(0 Example 5: angle between blades on 
conic flow surface Ad, 18°; other parameters 
same as standard example. 


CJ 


Relative Afaah number 
(ficfi fious T see sub legend) 



(g) Example 6: incompressible fluid; other parameters same as example 5 (Sg. 8 ( 0 ). Note 
that <p and Air are based upon constant c, T which if considered equal to speed of sound at 
inlet conditions of example 5 enables comparison of compressible (example 5) and incom- 
pressible (example 6) solutions for same impeller-tip speed, weight-flow rate, and so forth. 
For incompressible fluids the relative Mach number is fictitious and Is equal to relative 
velocity g divided by constant c«. 


Fiocri 8.— Concluded. Lines of constant Mach number relative to impeller. 


the 45° line on this plot. The error is, with very few excep- 
tions, less than 0.01, where the error is defined as 


Error=(y- 

w. 



The negative values of shown on the correlation plot 
correspond to eddies, which were shown in the previous 
section to form on the positive blade surface at low flow 
rates and at high tip speeds. The stream-function ratio is 
always zero at the positive blade surface and always 1.0 at 
the negative blade surface, 

Radial-velocity ratio V. — The radial-velocity ratio can be 
estimated by the following correlation equation (appen- 
dix B): 


y = V „ + ^1 ( F— F*) , + 

A(2 V - 1)1-^ TAfl- 


l 


[Ae- (AS) J -f (R—R,) {M t AS) , J (B 1 3 ) 


The estimated values of the radial-velocity ratio V ob- 
tained from equation (B13) are compared in figure 10 with 
the values of V obtained from the relaxation solutions of 
this report. Values are plotted for every other grid point 
across the passage at the values of f indicated by the sym- 
bols. Perfect correlation corresponds to the 45° line on this 
plot. The eiTor is less than 0.01 where the error is defined 
as 

Error = (V'—V) 


The negative values of V shown in the correlation plot 
correspond to the eddies that form on the positive blade 
surface at low flow rates and at high tip speeds. 

Tangential-velocity ratio U. — The relative tangential- 
velocity ratio can be estimated by the following correlation 
equation (appendix B): 


U'=AU S 


(Bl) 
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Figure 9.~Comparison between exact (relaxation solution) and estimated, (correlation equation (BIB)) values of stream-function ratio for various values of 72 (function of £ and ,19) and 0. 


The estimated values of the relative tangential-velocity 
ratio V obtained from equation (Bl) are compared in fig- 
ure 11 with the values of U obtained from the relaxation solu- 
tions of this report. Values are plotted for every grid point 
across the. first half of the passage for the values of £ indicated 
by the symbols. Perfect correlation corresponds to the 45° 
line on this plot. Except for the incompressible solution 
(example 6), the error is less than 0.01 where .the error is 
defined as ...... : ■ V ' 

Error = (U' — U) 

The relative tangential-velocity ratio, is always zero along 
the blade surfaces. 

Other flow conditions.— Other flow conditions within, the 


impeller can be, determined from the values of U and V 
obtained by the correlation equations (Bl) and (B13). 
The relative velocity ratio Q is given by 

Q2=U‘+V‘ (13) 

From the steady-flow energy equation, the temperature 
ratio is given by (equation (10), reference 1) 

" IWMry-Q 1 } 

where T is the absolute, static temperature of the gas and 
where the absolute whirl ahead of the impeller is assumed to 
be zero.. 
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Figure 10.— Comparison between exact (relaxation solution) and estimated (correlation equation (B13)) values of radial-velocity ratio for various values of R (function of { and ay) and 9. 


The density ratio p/p„ and the pressure ratio plp 0 are 
obtained from the temperature ratio by 


and 



i 

[(£ d / r ) 2 - e ]} 7 ” 1 


T 

[( f ? 3 / r ) 2 -< 2 2 ]} 7-1 


(12; 


The local relative Mach number is related to Q and the 
temperature ratio by 

M=Q^ 


Correlation of velocities along blade surfaces. — Of special 
interest, because of boundary-layer considerations, are the 
velocity distributions along the positive and negative sin- 
faces of the impeller blades, Q p and Q*. Estimated values 
of Qh obtained from the correlation equations are given by 
the curves in figure 12 for each of the numerical examples 
in this report. The values of Q K obtained by the relaxation 
solution are shown by the plotted points. The agreement is 
seen to be excellent in all cases. 

Estimated values of Q p obtained from the correlation 
equations are given by the curves in figure 13 for each of 
the numerical examples in this report. The values of Q p 
obtained by the relaxation solution are shown by the plotted 
points. Again the agreement is excellent. . 
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Figure 11.— Comparison between exact (relaxation solution]) and estimated (correlation equation (B 1)) values of relative tangential- velocity ratio for various values of R (function of £ and SO) 

and s. . 


Slip factor.— The impeller slip factor n is defined as the 
ratio of the absolute tangential velocity of the fluid at the 
impeller tip to the tip speed of the impeller. The following 
correlation equation is developed in appendix C for the slip 
factor of impellers with straight blades: - 

- (C2) 


For the standard solution, (A0), is 0.2095 radian and n, is 
0.934 so that the slip-factor equation becomes 

-" r - V = 1-O.315A0 (17) 

This equation is plotted in figure 14 together with Stodola's 
equation, which for straight blades is given by 

O.5OOA0 (IS) 
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The slip factors obtained for the numerical examples of 
this report are plotted as points in figure 14. These points 
indicate that the slip factor is independent of the flow rate, 
impeller-tip speed, and variation in passage-height ratio 
with radius and depends only on the angle between blades 
on the conic flow surface A 6. This conclusion is in agree- 
ment with Stodola’s equation although the magnitude of the 
slip factor given by Stodola is lower. 

APPLICATION OF CORRELATION EQUATIONS 

The correlation equations presented in the previous section 
are important because they provide rapid solutions (in the 
regions investigated herein, fig. 6) to the differential equa- 
tion, which determines the flow conditions (U, V, pfp e , and 
so forth) in any impeller (with straight blades) for all operat- 
ing conditions. An important application of the correlation 
equations is the determination of velocities along the flow 
surfaces because these velocities are significant in the study 
of the boundary layer. As examples of the information pro- 
vided by the correlation equations for such a study, the 
velocity ratios along the blade surfaces have been computed 
over a wide range of impeller-tip Mach number M r , flow 
coefficient <p, and angle between blades on the conic flow 
surface AS. These computations are presented in figures 15 
to 17. 

Impeller-tip Mach number M T - — The effect of M T on the 
velocity ratios along the positive and negative blade surfaces 
is shown in figure 15. In this figure all design and operating 
conditions (other than M T ) were maintained constant at the 
standard values. 

For M r equal to zero, the velocities are equal on both 
surfaces of the blade. For all other values of M r , the 
relative velocities are higher on the negative blade surface 
than on the positive blade surface and as the impeller-tip 
speed increases the difference in velocities along the two 
surfaces increases. Except near the tip, this increase in 
velocity difference results primarily from a decrease in 
velocity along the positive blade surface and for high values 
of M t this velocity becomes negative, which indicates the 
presence of an eddy. The small effect of M T on Q x results 
from a combination of effects. At higher values of M r , the 
difference between Q x and the mean radial-velocity ratio 
V m increases, but V„ itself decreases because of the increased 
density; the net result is only a small change in Q„ with 



Figure 14.— Comparison among slip factors obtained from relaxation solutions, correlation 
aquation UP, and Stodola’s equation (IS). 
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Figure 15.— Effect of impeller-tip Mach number on relative velocity ratio along positive 
and negative blade surfaces. Correlation equation (B13); flow coefficient p, 0.5; constant 
' flow area (m,— 1.0); angle between blades on conic flow surface AO, 12°; compressible 
flow (y, 1.4). ...... . . ... ... ... 



Figure 16.— Effect of flow coefficient on relative velocity ratio along positive and negative 
blade surfaces. Correlation equation ( B 1 3) ; impeller-tip Mach number Mr, 1.5; constant 
flow area (m,— 1.0); angle between blades on conic flow surface M, 12°; compressible 
flow (y, 1.4). ’ 


changes in M r . At the blade tip, the velocity becomes equal 
on both surfaces and this velocity decreases with increasing 
M T because of the. higher gas density. 

Flow coefficient <p . — The effect of flow coefficient ^ on the 
velocity ratios along the positive and negative blade surfaces 
is shown.in figure 16. In this figure, all design and operating 
conditions (other than the. flow coefficient) were maintained 
constant at the standard values. 

At each radius ratio, the difference between Q„ aDd Qj, is 
independent of f (that is, remains constant). The mean 
radial-velocity ratio V m , however, decreases with decreasing 
flow coefficient and for low flow coefficients Q p becomes 
negative, which indicates, the presence of an eddy. At the 
blade tip, the velocity becomes equal on both surfaces and 
this velocity increases with increasing flow coefficient because 
of the increased mean radial velocity V m . 

For ip equal to zero, the velocities arc equal on both 
surfaces of the blade (but opposite in sign) . As a result the 
pressures on' both blade surfaces are equal and no work is 
done by the impeller. The entire flow within the passage 
is an eddy. 

Angle' between blades on conic flow surface AS. — The 
effect of AS on the velocity ratios along the positive and 
negative blade surfaces is shown in figure 17. In this figure 
all design and operating conditions other than A9 were 
maintained constant at the standard values. Increased 
values of AS indicate fewer blades. 

The correlation equation (B13) used to obtain the curves 
plotted in figure 17 is developed in terms of the transformed 



Figure 17.— Effect of angle between blades on conic flow surface AO on relative velocity 
ratio along positive and negative blade surfaces. Correlation equation (B13); dashed 
lines obtained from simplified analysis (equation (10)); flow coefficient *>, 0.5; Impeller-tip 
Mach number Afr. 1.5; constant flow area (m, -1.0); compressible flow (7, 1.4). 
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coordinate £, which is related to the radius ratio R by an 
expression (equation (B3)) that includes the angle between 
blades on the conic flow surface A0. The minimum value 
of £ for wliich standard values of U, V, and so forth are 
given in tables I to III is —1.7030, wliich for values of AS 
less than standard (AS= 12°) corresponds to values of R 
greater than 0.7, as indicated by equation (B3) and shown 
in figure 17. For values of R less than those resulting from 
£ equal to — 1 .7030, however, the simplified analysis presented 
in reference 1 for impellers with straight- blades may be used 
to extrapolate the curves to all lesser values of R. In 
terms of the mean velocity ratio V m (equation (B9)), the 
simplified equation for the velocity ratio becomes (Note 
that the simplified analysis of reference 1 assumes U— 0.) 

Q=r=r„+flA/rAfl[2(£)-l] (19) 

where V equals V m midway between blades '(that is, when 
0=A0/2). Equation (19) has been used to extrapolate the 
curves for values of A0 less than 12° in figure 17 (dashed 
lines). For values of £ less than —1.7030, the simplified 
analysis given in reference 1 can be used to determine the 
flow conditions within the impeller passage for any design 
and operating condition. For values of £ greater than 
— 1.7030, the methods of this report must be used to deter- 
mine V, pjp„, and so forth. 

For A0 equal to zero, the velocities are the same on both 
surfaces of the blade and are equal to the mean velocity 
V m (dashed line). This mean velocity is the same for all 
values of A0, but the difference between I r m and Q p increases 
with increasing values of A0 so that, for large values of A0, Q v 
becomes negative, which indicates the presence of an eddy. 

SUMMARY OF RESULTS AND CONCLUSIONS 

Six numerical examples are presented for steady, two- 
dimensional, compressible, nonviscous flow in centrifugal 
compressors with thin straight blades, the center lines of 
which generate the surface of a right circular cone when 
rotated about the axis of the compressor. A seventh ex- 
ample is presented for incompressible flow. The solutions 
were obtained in a region of the compressors (including the 
impeller tip) that was considered to be unaffected by the 
inlet configuration of the impeller and by the diffuser vanes. 
(That is, the impeller inlet and the diffuser vanes, if any, 
must be far enough removed from the region investigated 
not to affect the flow appreciably in that region.) The 
effects of variations in the following parameters were investi- 
gated: (1) flow coefficient (flow rate), (2) impeller-tip Mach 
number, (3) exponent for the variation of passage-height 
ratio with radius ratio, and (4) the angle between blades on 
the conic flow surface. Each solution applies to radial- and 
mixed-flow compressors (and turbines with the rotation and 
flow direction reversed) with various cone angles but with 
the same angle between blades on the conic flow surface. 
The numerical results are presented in plots of the stream- 
lines and constant Mach number lines. 

Correlation equations are developed whereby flow condi- 
tions (streamlines, velocities, pressures, impeller slip factor, 
and so forth) within any impeller with straight blades can 
be determined (for flow region investigated) from the flow 

*i -i , > 


conditions of the standard solution presented. As examples 
of the information provided by the correlation equations, the 
velocities along the blade surfaces have been computed (and 
plotted) over a wide range of impeller-tip Mach number, 
flow coefficient, and angle between blades on the conic flow 
surface. 

The principal conclusions resulting from the work pre- 
sented herein are: 

1. The exponent m, which was used in this analysis to 
specify the variation in passage-height ratio with radius 
ratio, has only a small effect (for practical values of m) on 
the area variation in the flow region investigated and there- 
fore has only a small effect on the flow in this region. 

2. An eddy forms on the positive blade surface at high 
impeller-tip Mach numbers, low flow coefficients, and large 
angles between blades on the conic flow surface. At an 
impeller-tip Mach number of 2.0 (and the standard values 
of flow coefficient and angle between blades on the conic 
flow surface), the eddy occupies more than 50 percent of 
the flow area at a radius ratio of 0.90. 

3. Compressibility has a large effect on the streamline 
configuration within the compressor (and therefore on the 
other flow conditions). For example, the large eddy that 
exists in the compressible-flow example (with the same param- 
eters as the incompressible example; completely disappears 
in the incompressible example. 

4. The maximum- relative Mach number occurs on the 
negative blade surface at a radius ratio well within the im- 
peller, and the flow- decelerates along the surface of the blade 
from this point to the blade tip. This deceleration, which 
for impellers with straight blades and with the usual type of 
area variation with radius ratio becomes rapid near the 
blade tip, is conducive to boundary-layer separation. 

5. If the boundary-layer wake in the vaneless diffuser is 
neglected, the flow conditions in the vaneless diffuser follow- 
ing the impeller become essentially uniform at a value of £ 
approximately equal to 0.45. 

6. For the high impeller-tip Mach numbers investigated 
(and if boundary-layer effects are neglected), the relative 
velocities at the impeller tip are low because the high 
impeller-tip Mach numbers result in high fluid densities. 

7. The maximum relative Mach number (on the negative 
blade surface) is increased by increasing the flow coefficient 
or the angle between blades on the conic flow surface (that 
is, the number of blades), but it is little affected by the 
impeller-tip Mach number. 

8. The impeller slip factor is independent of the impeller- 
tip Mach number, compressor flow coefficient, variation in 
passage-height ratio with radius ratio, and compressibility of 
the fluid. The slip factor is a function only of the angle 
between blades on the conic flow surface. 

9. The difference between velocities on the positive and 
negative blade surfaces increases with increasing tip Mach 
number and angle between blades on the conic flow surface 
but is independent of the flow coefficient. 

Lewis Flight Propulsion Laboratory, 

National Advisory Committee for Aeronautics, 
Cleveland, Ohio, June 21, 1949. 



APPENDIX A 

SYMBOLS 


The following symbols arejised in this report: 

A ratio (equation (B2)} 

a annular flow area (normal to conic. flow sur- 
face) ' ;■ : 

B number of passages (orblades) 

b grid spacing (fig. 3) 

e ..local speed of sound 

exp exponential, [exp (x)=e I ) 

II passage-height ratio, hjh T .......... 

h passage height (normal to conic flow surface) 

M relative Mach number. _ 

Mr impeller-tip Mach number (equation (9)) 

m passage-height exponent (equation (11)) 

p absolute static pressure . . 

Q relative velocity ratio, qjc 0 

q velocity of fluid relative to blades, ti 2 +r ! 

B conic-radius, ratio (coordinate of conic, flojv 

surface) , rjr T 

R residual ........ - : 

r conic radius - (distance along conic element 

from apex of cone) (fig. 1) 

T absolute static (stream) temperature' 

U relative tangential-velocity ratio, u/c 0 

u tangential component of* q (in positive direc- 

tion of 8 ) (fig. 2) 

V . radial-velocity ratio, vlc 0 

v radial component of q (along conic element) 

(fig. 2) . ■ 

W total compressor flow rate 

a cone angle (fig. 1) \ 


y 

v 

8 


A8 


A* 

P 

<P 


Subscripts: 

m 

n_ ' ~ 

0 “ 

V . 

8 ’ , 

T 

1, 2, 3, and 4 
Superscripts: 


ratio of specific heats 
transformed coordinate (equation (10b)) 
angle (coordinate of conic flow surface, posi- 
tive in counterclockwise direction), radians 
unless otherwise specified 
angle between blades on conic flow surface 
(equation (6), (fig. 6)), radians unless other- 
wise specified 
impeller slip factor 

. transformed coordinate (equation (10a)) 
weight density of fluid 

compressor flow coefficient (equation (7)) 

dimensionless compressible stream function 
(equations (4a) and (4b)) ... _ 

impeller angular velocity (in positive direc- 
tion of 8) - ■ 

mean value at given radius ratio 
negative blade surface (blade surface in nega- 
" tive direction of 8) 

absolute stagnation condition in region of 
uniform flow upstream of impeller 
positive blade surface (blade surface in posi- 
tive direction of 8) 

standard solution . , ...... , . 

impeller tip 

grid points adjacent to point being considered 
(fig- 3) 

estimated value 


APPENDIX B , 

CORRELATION EQUATIONS 


Correlation equations are developed whereby flow con- 
ditions (U, V, and within any impeller with straight 
blades can be determined from the flow conditions for the 
standard solution of this report. 

TANGENTIAL-VELOCITY RATIO U 

A plot of U against 17 (equal to 8/A8, equation (10b)) at $ 
equal to 0 for various design and operating conditions (used 
in the relaxation solutions of this report) is shown in figure 18. 
These curves are representative examples of the variation in 
U for all values of £ at which U is significant. For a given 
value of £ and 17, tire tangential-velocity ratio U is seen to 
be a function primarily of the impeller-tip Mach number 
M t and the pngle between blades on the conic flow surface AS, 
The tangential-velocity ratio is essentially independent of 
the compressor flow rate (flow coefficient <p) and the vari- 
ation of passage-height ratio with R (passage-height ex- 
ponent m) . This dependence of U on M T and A 8 only (for 
a given value of £ and y) was found to exist at all radius 
ratios at which U is significant and was found to be a direct 

relation such that . . ... : . . 

U=KM t A8 


where K is a function of £ and y that is constant for all 
design and operating conditions. In terms of the standard 
solution, therefore, 


(Bl) 

(B2) 

and where the subscript s refers to values from the standard 
example and the prime indicates the estimated value for the 
nonstandard example. Equation (Bl) is the correlation 
oquationfor the relative tangential-velocity ratio. 

The correlation equation (Bl) and the other correlation 
equations to be developed in this appendix refer to the same 
transformed coordinates £ and y for both the standard and non- 
standard. quantities in the equations. From equations (10a) 
and (10b), these transformed coordinates are related to the 
coordinates in the physical (R,8) plane by 


where 


U'=AU, 

, M t A8 
(M t A8), 


and 


R=exp (£Afl) 


8=yA8 


m 

(B4) 
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where V m is the mean radial-velocity ratio, which is obtained 
from continuity considerations as follows: 



Fiook* 18. — Variation in relative tangential-velocity ratio across impeller passage for £— 0. 
Plotted points obtained from relaxation solution. 

RADIAL-VELOCITY RATIO V 

The correlation equation for the radial-velocity ratio is 
obtained from the equation for irrotational absolute motion 
of a fluid particle. From reference 1, 

■nr -U,*uibV 

-^T-R-t b R R^g 


(B5) 


which, after multiplying by A0, becomes 

Plots of V against ij for various design and operating 
conditions and for £ equal to 0 and —1.0 are shown in fig- 
ure 19. These plots are representative of the variation in \ 
at all values of £ within the impeller. The slopes of the 

velocity profiles are seen to be nearly constant (except in 
the immediate vicinity of the blade tip) so that from equa- 
tion (B5) 

(B6) 


(B7) 


and equation (B5) can be integrated to give 
T = T j-)-2i?il/r(A0)ij-t- R{A6) ^ t" 

where V equals V, and ij equals zero along the positive blade 
surface. 

The velocity ratio V „ can be evaluated from the condition 
V=V m 

when 


'-5T 0 - 6 ' 


(B8) 


W=p m V m rhBA8 


or 


W 


<p 


Po&tCo Pe 


VJtH 


so that 


^RH 

Po 


m 


The mean density ratio pmlPo in equation (B9) is obtained 
from equation (12) by assuming that for straight blades the 
relative tangential velocity is zero and the radial-velocity 
ratio is equal to the mean radial-velocity ratio V m 







Fiqcbi 19.— Variation in radial-velocity ratio across Impeller passagd for two values of f. 
Plotted points’obtained from relaxation solution. 
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Equation (12) becomes 
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Hi- 

Po ) 


7—1 


(RMt) 




(BIO) 


Equation (BlO) gives an average value of the density ratio, 
which is assumed to be a function of £ (that is, (log s E)/A0) 

and independent of t; (that is, 6/A6). 

From equation (B7) and the condition given by equa- 
tion (B8) , the velocity ratio V v becomes. 

T T , . , , RA6/U . bU\ 

v„=v m -RM T A8 — 2 “ {ji+M ) ::: 

and equation (B6) becomes 

V=V m +RM T A8(2r,~l) +~4 6 (jl + §f) (2u-l) 

From equation (B3), ■ ; - 

bU = bUdl-__ 1 bV 
bR £>£37? RB b£ 


(Bll) 

V 


and from equation (Bl), 

f'U,bU\ A ( TT Afl , bU,\ ' 

\R + bRrm{ u ' Ad+ ~w) 

so that, equation (Bll) becomes 

V=V m +RM r A0(2„-l)+| (W b j)(2^-1) 
which is solved for b UJbt; 


(B12) 


bU, 


V— V 


2RM r A0 




i (2u-l) 


• U,A6 


The term bU s jb% can be eliminated by equating it for the 
standard and nonstandard cases. Therefore, 

v'=v m +A(v~v m ),+ . . , 

-4(2ij— 1) [A0-(Atf) t ] + (R—R,) (M T Afl) ,| (B13) 

Equation (B13) is the correlation equation for the radial- 
velocity ratio. ■;■■■■■ ' ;; ; - 


STREAM-FUN CTION RATIO 

The stream-function ratio is obtained from the definition 
of the stream function and from the correlation equation for 
the radial-velocity ratio T 7 . The definition of the stream 
function is.given.by equation (4a): 


oe p 0 


(4a) 


which, after being multiplied by AS and combined with 
equation (.Bll), becomes . 

^^RH(A^V m +RM r {A0){2 v -l)+^(^+^{2n-l)^ 




(B14) 


.From the_condition given by equation (B6) and with the 
density ratio replaced by its mean value for £ (that is, 
(log„R)/A£?), equation (B14) can be integrated to give 

r.-k [ r - +BMA 4 ® ^ + ~r (k + sl) '’- 1 *] ’ 

(BIS) 

where ipn is given by equation (5) and is equal to zero 
when ij is equal to zero. Introducing equation (B12) and 
solving for & £/,/d£ give 


bU, 

■sr 


V „ 


q-o- 


RM r (A9) ( 17 — 1) 17 


^( 17 — 1)17 


■U,AQ 


The term bUJbi, can be eliminated by equating it for the 
standard and nonstandard cases; therefore - 

jl) [aS-{i9),]-|- (H-BJ (J/rAOj. j) (BIS) 

Equation (B16) is the correlation equation for the stream- 
junction ratio. 


APPENDIX C 

CORRELATION EQUATION FOR IMPELLER SLIP FACTOR 


The impeller slip factor is defined as the ratio of the average 
absolute tangential velocity of the fluid leaving the impeller 
tip to the tip speed of the. impeller. This definition results 
in the following equation (reference 1): 


/»=! + 


U m 

M v 


(Cl) 


Also, from equations (Bl) and (B2) 

TL± = (TJ m \ A 9, 
Air \MrJt (AS) , 

which, from equation (Cl), becomes 


or 


A — i = (u,— 1 ) 

ti' — l — (1— fi t ) 


a e 
(A B)~ 

a e 
(A e). 


(C2) 


Equation (C2) is the correlation equation for the slip factor. 
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TABLE I— STREAM-FUNCTION RATIO FOR STANDARD SOLUTION 


[Conditions for standard example: (®),=0.5; CMr),= 1.5; (n0,=— 1.0; (i$),=0.20941 radian; Ct),=IAI 


0»), 


i 

* n (equals 8[A9) 


D 

0.1 

0.2 

0.3 

... 

0.5 

0.6 

0.7 

0.8 

0.8 

1.0 

0,3163 

a to 

-1.7030 

n 

a oio 

0.095 

0.165 

0. 244 

0.341 

0.451 

0.570 

0.702 

0.843 

1.000 

. 2041 


-1. 3736 


.031 

.078 

.143 

.220 

.316 

.427 

.519 

.686 

.810 

1.000 

.2735 


- L 0654 


.022 

.061 

. 120 

.187 

.291 

.401 

.529 

. 871 

.831 

1.000 

.2530 

,85 

-.7760 


.014 

.047 

.101 


.270 

.383 

- .509 

.655 

.824 

3.000 

. 23 IS 

.90 

-..5031 


.011 

.042 

.067 


.255 

.377 

.503 

.650 

.819 

1.000 

.2275 

.92 

— . 39 S 1 


.012 

.016 

. 102 


.270 

.382 

.506 

.651 

.816 

1.000 

.2202 

.94 

-. 2954 

0 

.016 

.054 

.113 

.189 

.283 

.392 

.517 

.660 

.820 

L 0 OO 

.2133 

.90 

-.1019 

0 

.023 

.069 

. 132 

.212 

.305 

.413 

.536 

.675 

.828 

LOGO 

.2100 

.97 

-.1154 

■■ 

.029 


.146 

.228 

.321 

.429 

.550 

.684 

.832 

L 000 

.2067 

.98 

-.0965 

■9 

.038 


.163 

.247 

.341 

.447 

.566 

.697 

.839 

1.000 

. 2034 

.99 

-.0130 

■■ 

.048 

ttf 

.185 

.270 

.366 

.472 

.586 

.713 

.850 

L 000 

.2002 

1.00 

0 

u 

.062 

HI 

» 210 

.297 

.393 

.499 

.611 

.734 

.862 

■ 1.000 


TABLE II — RADLA.L- VELOCITY RATIO V FOR STANDARD SOLUTION 


[Conditions for standard example: M,*0.5; LVr).“I.S; Cm),*— L0; CA9), =0.20911 radian; C7),=t.4] 


G »*, 

R . 

t 

n [equals 9/ A 9) 

0 

0.1 

0.2 

0.3 

0.4 

0.5 

0.6 

0.7 

0.8 

0.9 

1.0 

0.3163 

0.70 

- I . 7030 

0.103 

0- 140 


0.234 

0.278 

0.323 

0. 368 

0.412 

0.456 

0.501 

0.549 

.2941 

. 75 

-1.3736 

.065 

.112 


.206 

.253 

.301 

.346 

.396 

.443 

.492 

.513 

.2735 

.SO 

— L 0654 

.033 

.082 


.181 

230 

.280 

.329 

.378 

.428 

.479 

.533 

.2536 

.&5 

-.7760 

.008 

.059 

.110 

.160 

.210 

,259 

.309 

.357 

.108 

.460 

.515 

. 23 I S 

.90 

-.5031 

-.001 

.049 

.099 

.148 

.194 

.239 

.286 

.330 

.378 

.430 

.483 

.2275 

.92 

-.3981 

.002 

.051 

.100 

.146 

.189 

.231 

.273 

.315 

.363 

.413 


.2202 

.94 

29-54 

.010 

.058 

.105 

.148 . 

.186 

.224 

.261 

.299 

.342 

.390 

.433 

.2133 

.96 

-.1949 

.025 

.072 

.115 

.151 

. 184 

.210 

.248 

.280 

.316 

.360 

.407 

-2100 

.97 

-.1454 

.037 

.082 

.122 

.155 

. J 84 

.213 

.241 

.270 

.301 

.342 

.388 

.2067 

.98 

-.0965 

.055 

.095 

.131 

. I 5 S 

.184 

.209 

.234 

.259 

.286 

.321 

.362 

,2034 

.99 

— - O 4 S 0 

.091 

.115 

.139 

.161 

. 1 S 4 

.206 

.227 

.249 

.271 

.293 

.316 

.2002 

1.00 

0 

.199 

.131 

.118 

.166 

. 1 S 4 

.202 

.220 

.238 

.254 

.271 

.199 

‘ 


TABLE III— RELATIVE TANGENTIAL-VELOCITY RATIO U FOR STANDARD SOLUTION 

[Conditions for standard example ; f . UrL — JLS ; (A5),=0.20<H4radianI 


































































